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We establish a connection between the distribution functions of quarks in spin-one hadrons and 
the helicity matrix for forward scattering of antiquarks off spin-one hadrons. From positivity of 
this matrix we obtain inequality relations among the distribution functions. Analogous relations 
hold also for fragmentation functions. The bounds we obtained can be used to constrain estimates 
of unknown functions, occurring in particular in semi-inclusive deep inelastic scattering or e + e~ 
annihilation with vector mesons in the final state. 
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In recent years some attention has been given to distribution functions characterizing spin-one targets, starting from 
the work of Hoodbhoy, Jaffe and Manohar [jjj. Unfortunately, the only available spin-one target is the deuteron, which 
is essentially a weakly bound system of two spin-half hadrons. In the approximation of independent scattering off the 
nucleons, the specific spin-one contribution to the deuteron structure functions is expected to be small, except maybe 
for the low-x region |^|§. New precision measurements of the deuteron structure functions coming from HERMES 
\^ ' may provide an experimental test of this expectation, as already observed in several papers [§-(6) . 

. Instead of pointing the attention to spin-one targets, it is also possible to analyze spin-one final state hadrons 
t-H ' in e + e~ annihilation or in semi-inclusive deep inelastic scattering. This idea was first considered by Efremov and 
Teryaev 0. A systematic study was accomplished by Ji |8|, who singled out two new fragmentation functions, the 
function b\, in analogy to the distribution function b%, and the time- reversal odd (T-odd) function h-. These new 
fragmentation functions can be observed in the production of vector mesons, e.g. p, K*, <f>. However, these functions 
require polarimetry on the final-state meson, which can be done by studying the angular distribution of its decay 
I products (e.g. tt + tt~ in the case of p° meson). In this sense, vector meson fragmentation functions represent just a 
specific contribution to the more general analysis of two-particle production |9|,|l0|| near the vector meson mass. 
O j| In a recent work [Tl| ] we contributed to the study of spin-one distribution and fragmentation functions by defining all 
D . possible functions occurring at leading order in 1/Q when also partonic transverse momentum is included, extending 
1 the work done for spin-half hadrons in jl^] . 

The study of spin-one fragmentation functions has a twofold relevance: it offers new insights in the physics of 
mesons and it serves as a probe of particular nucleonic distribution functions. With regards to the second issue, we 
point out that for 1-particle inclusive leptoproduction of a spin-one meson, at leading (zeroth) order in 1/Q there are 
■ five chiral-odd fragmentation functions that can be observed in combination with the quark transversity distribution 
function ]l3j , the measurement of which is at present attracting the attention of experiments such as HERMES, 
COMPASS and RHIC. 

In this situation, it seems to be extremely useful to gather as much information as possible about spin-one dis- 
tribution and fragmentation functions. As a first step in this direction, we present here positivity bounds on them, 
following what has been done in Ref. |l4|] for the spin-half case. 

II. BOUNDS ON TRANSVERSE MOMENTUM INTEGRATED FUNCTIONS 

Distribution functions appear parametrization of the light-cone correlation function |p^-|r7|| 

(1) 



J+=Jt=0 



* Corresponding author. E-mail:bacchett@nat .vu.nl 



1 



depending on the light-cone fraction x = p + /P + , where P is the momentum of the target hadron andp the momentum 
of the outgoing quark. Since we are dealing with spin-one hadrons the polarization state of the target must be specified 
using a spin vector S and a spin tensor T . In the target rest-frame, we can conveniently parametrize the spin vector 
and tensor as [] 
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At leading twist (leading order in l/Q) in inclusive leptoproduction, the relevant part of the correlator is $7^ 
transpose represents the forward amplitude for antiquark-hadron scattering, M = ($7 + ) . 
scattering amplitude, the matrix M in the parton <S> hadron spin space is positive definite, 
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where ip + = V+ijj = ^7 7 + V' is the good component of the quark field |Q and where A, A' = 1,0, 
For spin-one hadrons, $(a;)7 + is parametrized in terms of five distribution functions 

Hxh+ = {/!(*) + gi {x) S L 75 + hi(x) 75)&t + Hll{x) S ll + i h 1LT (x) p LT }v+. 
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The last distribution function is constrained to be zero by time-reversal invariance. Nevertheless, we choose to keep it 
because in the analysis of fragmentation functions the analogous contribution, Hxlt(z), cannot be constrained in the 
same way. To follow a more systematic naming of the functions when we will include transverse momenta, we felt the 
need to change the names of the transverse momentum independent functions. We are aware that changing notations 
often poses some undesired difficulties, but we believe in the convenience of using a notation that harmoniously 
connects transverse momentum dependent to transverse momentum integrated functions. With the parametrizations 
in Eq. ([}]) and Eq. ([5]), however, our function Jill is precisely the same as the function 61 of Hoodbhoy, Jaffe and 
Manohar 0. 

To construct the matrix M for the leading order correlation function, only two basis states in Dirac space are 
relevant, corresponding to good components of the right and the left handed partons. This is particularly transparent 
when writing Eq. (Q) in chiral representation. Therefore, we can effectively reduce the four-dimensional Dirac space 
to the 2x2 good parton chirality space. In other words, only the scattering of good antiquarks off hadrons gives a 
non vanishing contribution to leading twist. 

The leading-twist scattering matrix takes then the form 
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The result contains the dependence on the hadron spin vector and tensor. We aim at obtaining a 6 x 6 matrix 
encompassing also the hadron spin space, with no dependence on the spin vector and tensor. To this goal, we need 
to use the hadron spin density matrix, p(S,T), and invert the relation 



x Note that the the parametrization of T is different from the one of Ref. There is a factor —2/3 difference between the 
old definition of Sll and the new, so that now —1/3 < Sll < 2/3. 
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M ZJ (x;S,T) = pyx(S,T) M iX , jX ,{x). 
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The spin density matrix for spin-one hadrons can be decomposed on a Cartesian basis of 3 x 3 matrices, using the 
(rank-one) spin vector, S l , and the rank-two spin tensor, T lJ , 
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where the matrices X 1 are the generalization of the Pauli matrices to the three-dimensional case and where we chose 



S v = I ( E <j)J + - \l 6 ij . 

Inverting Eq. (0), we can eventually reconstruct the 6x6 leading- twist scattering matrix 
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where all the functions on the right-hand side depend only on x. 
From the positivity of the diagonal elements we obtain the bounds: 

fi(x)>0 
-\h{x)<h LL {x) < 3/x(x) 

\gi(x)\<h(x)-±f 1LL (x) < \h{x), 

while positivity of 2-dimensional minors gives the bound p9| 

Ihix)} 2 + [h 1LT {x)f <\[h{x) + lf ltL {x)\ h{x)+ 9l {x)~ 
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This bound is a generalization of the Soffer bound and must be fulfilled by any spin-1 target. When /iilt = due to 
time-reversal invariance, the bound can be reduced to 
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An analogous calculation can be performed for fragmentation functions, which describe the hadronization process 
of a parton into the final detected hadron. In this case the transverse momentum independent correlator is |20|] 



A ^) = E / ^ 7 ^e ik <(0\MO\Ph,S hl T h ;X)(P h ,S h ,T h ;X\i, J (0)\0) 
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depending on the light-cone momentum fraction z = P h jk , where Ph is the momentum of the outgoing hadron and 
k the momentum of the fragmenting quark. 
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The analysis of the leading-order scattering matrix can be developed in complete analogy to the distribution 
correlation function. The matrix M — now represents a hadron-quark decay matrix, which is again positive 

definite. Instead of using the hadron spin density matrix, we should employ the hadron decay matrix, expressing it 
in terms of vector and tensor analyzing powers Al, At, All, etc. The final result, however, would still have the form 
of Eq. (|Tc|), with the z-dependent fragmentation functions D\,G\,D\ll,H\ and H\lt replacing the x-dependent 
distribution functions fi, gi, /ill, hi and Hilt- Despite the differences in notation, our fragmentation functions 
correspond to the ones of Ji II (i.e. Dill = &i and Hilt — hj). Note that in the context of a fragmentation process 
it is not justified to discard T-odd functions. We point out that the description in terms of the matrix in spin space 
allows a connection to the analysis performed in . 

It is particularly useful to set bounds on the function Hilt, which is a possible candidate to probe the transversity 
distribution of the nucleon in inclusive deep-inelastic vector meson production. It appears in connection with the 
transversity distribution in a (sin (tp^n + 4>s)) transverse spin asymmetry, where 4>s and are the azimuthal angles 
with respect to the lepton scattering plane of, respectively, the spin of the target and the difference of the momenta 
of the two decay particles |L3| . The equivalent of Eq. ([TI]) for fragmentation functions implies the bounds 
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D^z) + -D 1LL {z) D 1 (z)--D 1LL (z) < - [D^z)] 
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The less restrictive bounds are particularly relevant because in the specific case vector meson decay the vector analyzing 
powers are zero, making it impossible to measure the fragmentation functions H\ and G\ . Note that some information 
on the function Dill can be already extracted from experimental measurements [22 23] . 



III. BOUNDS ON TRANSVERSE MOMENTUM DEPENDENT FUNCTIONS 



To include the dependence on the parton transverse momentum, p T , the correlation function can be defined as |12|] 
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For convenience, the correlation function can be decomposed in several terms in relation to the polarization state of 
the target, i.e. $ = $[/ + $l + + &ll + $lt + &tt- To leading order in l/Q, these terms can be decomposed as 
(we identify T-odd terms by enclosing them in round parentheses) 
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Following steps analogous to the previous section, we can reconstruct the complete 6x6 scattering matrix. The 
inclusion of p T dependence makes all the entries of the matrix to be non-zero. It is also convenient to define the 
functions 
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In the rest of this Section, unless otherwise specified, all the functions are understood to depend on the variables x 
and pip. 

In principle, if we fully exploit the condition of the scattering matrix to be positive definite, we can write several 
relations involving an increasing number of different functions. We feel this to be an excessive task if compared to 
the exiguity of information we have on the involved functions. Therefore, here we choose to focus on the relations 
stemming from positivity of the two-dimensional minors of the matrix. 

Because of the symmetry properties of the matrix, only nine independent inequality relations between the different 
functions are produced. To simplify the discussion, it is useful to identify some of the T-odd functions as imaginary 
parts of the T-even functions, which become then complex scalar functions. The following replacements are required 
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Assuming the above relations, we can write the following inequalities: 
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A further simplification can be performed by introducing the positive real functions 

A{x,p 2 T ) = h{x,p 2 T )+ gi {x,p 2 T ) - ^fiLL(x,p 2 T ), 
2 

B(x,p T ) = fl(x,p T ) + - flLL{x,p T ), 

C{x,p 2 T ) = fi{x,p 2 T ) - gi{x,p 2 T ) - ^Hll{x,p%) 
and express the combinations of distribution functions occurring in the various entries of the matrix as 

hx = a^J\AB, 

9it + flLT = b V AB , 
9it - Jilt = c V£~B, 
hi T = dVCB, 
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where the coefficients a, 6, c, d, h, j are complex functions of the variables x and p\, and the coefficients e, /, g are 
real functions of the same variables. Time-reversal invariance in the distribution sector constrains the coefficients a, 
b, c, d, h, j to be real functions of the variables x and p\ and the coefficients e, /, g to vanish. 
The scattering matrix can now be expressed in a concise way as 
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Positivity of the two-dimensional minors of this matrix requires the absolute value of the coefficients a to j to be 
smaller or equal to unity. It is possible to consider higher dimensional minors to extract other useful bounds. The 
choice of the appropriate minors will depend on the availability of measurements of specific functions. 



IV. CONCLUSIONS 

In this letter, we studied the leading-order part of the correlation functions for parton distribution inside spin-one 
targets or parton fragmentation into spin-one hadrons. We cast the correlation functions in the form of scattering 
matrices in the parton chirality space ® hadron spin space. Positive definitcness of these matrices allowed us to set 
bounds to the involved distribution and fragmentation functions. These elementary bounds can act as important 
guidances to estimate the magnitude of otherwise unknown functions and, consequently, the magnitude of spin and 
azimuthal asymmetries appearing in several processes containing spin-one hadrons. In particular, these bounds are 
significant for one-particle-inclusive deep inelastic scattering and e + e~ annihilation with vector mesons in the final 
state. Among our results, we presented a bound that generalizes the Soffer inequality to the case of spin-one targets. 
Furthermore, we proposed a bound on the fragmentation function Hilt- Since the latter can appear in connection 
with the transversity distribution in a (sin {</>^ n + 4>%)) asymmetry, our bound can be useful in the extraction of the 
transversity distribution, where, for instance, it could be important to constrain the absolute normalization of the 
function. 



G 



ACKNOWLEDGMENTS 



This work is supported by the Foundation for Fundamental Research on Matter (FOM) and the Dutch Organization 
for Scientific Research (NWO). 



P. Hoodbhoy, R.L. Jaffe and A. Manohar, Nucl. Phys. B 312 (1988) 571. 

N. N. Nikolaev and W. Schafer, Phys. Lett. B 398, 245 (1997) [Erratum-ibid. B 407, 453 (1997)]; 

K. Bora and R. L. Jaffe, Phys. Rev. D 57, 6906 (1998). 

M. Anselmi no, A. Efremov an d E. Leader, Phys. Rep. 261 (95) 1-124. 



S. Kumario, tiep-ph/9302320 



A.V. Efremov and O.V. Teryaev, iep-ph/99 10555 



A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 36 (1982) 557 [Yad. Fiz. 36 (1982) 950]. 
X. Ji, Phys. Rev. D 49 (1994) 114. 

A. Bianconi, S. Boffi, R. Jakob and M. Radici, Phys. Rev. D 62 (2000) 034008. 

R.L. Jaffe, X. Jin and J. Tang, Phys. Rev. Lett. 80 (1998) 116; Phys. Rev. D 57 (1998) 5920. 

A. Bacchetta and P.J. Mulders, Pys. Rev. D 62 (2000) 114004. 

P.J. Mulders and R.D. Tangerman, Nucl. Phys. B 461 (1996) 197; Nucl. Phys. B 484 (1997) 538 (E). 

A. Bacchetta and P.J. Mulders, in preparation; Proceedings of R1KEN-BNL Research Center Workshop, Volume 29, 
BNL-52612. 

A. Bacchetta, M. Boglione, A. Henneman and P.J. Mulders, Phys. Rev. Lett. 85 (2000) 712. 

D.E. Soper, Phys. Rev. D 15 (1977) 1141; Phys. Rev. Lett. 43 (1979) 1847. 

R.L. Jaffe, Nucl. Phys. B 229 (1983) 205. 

A.V. Manohar, Phys. Rev. Lett. 65 (1990) 2511. 

J.B. Kogut and D.E. Soper, Phys. Rev. D 1 (1970) 2901. 

J. Soffer, Phys. Rev. Lett. 74 (1995) 1292. 

J.C. Collins and D.E. Soper, Nucl. Phys. B 194 (1982) 445. 

M. Anselmino, M. Boglione, J. Hansson and F. Murgia, Phys. Rev. D 54 (1996) 828. 
K. Ackerstaff et al. [OPAL Collaboration], Phys. Lett. B 412 (1997) 210. 
P. Abreu et al. [DELPHI Collaboration], Phys. Lett. B 406 (1997) 271. 



7 



